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Abstract:  The  deformation  of  boat  hull  bottom  panels  during  the  initial  phase  of  slamming  is 
studied  analytically  using  a  linear  elastic  Euler-Bemoulli  beam  as  a  representation  of  the  cross 
section  of  a  bottom  panel.  The  slamming  pressure  is  modeled  as  a  high-intensity  peak  followed  by 
a  lower  constant  pressure,  traveling  at  constant  speed  along  the  beam.  The  problem  is  solved 
using  a  Fourier  sine  integral  transformation  in  space  and  a  Laplace-Carson  integral  transformation 
in  time.  Deflection  and  bending  moment  as  functions  of  time  and  position  for  different  speeds, 
bending  stiffnesses,  etc.  are  given.  In  particular  the  effect  of  slamming  load  traveling  speed  on 
structural  response  of  the  simplified  bottom  structure  is  investigated.  It  is  found  that  rather  large 
deflections  and  bending  moments  are  encountered  at  certain  speeds  of  the  pressure,  which 
suggests  that  bottom  panels  may  benefit  from  tailoring  their  stiffness  and  mass  properties  such 
that  loads  are  reduced.  This  would  vary  with  boat  particulars  and  operation  (deadrise  angle,  mass, 
speed,  sea  state,  etc).  The  importance  of  the  high-intensity  pressure  peak  often  encountered  during 
slamming  is  also  studied.  It  is  seen  that  for  relatively  slow  moving  slamming  loads  the  pressure 
peak  has  little  influence.  However,  for  faster  moving  loads  its  influence  can  be  significant. 
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1 .  Introduction 


Some  of  the  highest  loads  on  high-speed  boats  are  due  to  bottom  slamming.  Slamming  pressures 
are  very  dynamic  and  vary  significantly  over  the  bottom.  Typically  slamming  starts  with  a  high- 
intensity  pressure  peak  that  travels  rapidly  over  the  bottom  from  the  keel  towards  the  chines.  The 
pressure  peak  is  usually  followed  by  a  lower  and  essentially  constant  pressure.  The  pressure  peak 
magnitude  and  propagation  speed  depend  heavily  on  the  impact  velocity  and  deadrise  angle  of  the 
boat.  Slamming  peak  pressures  have  been  experimentally  measured  to  reach  8  MPa  or  even  more 
(e.g.,  Faltinsen  [1]),  which  is  close  to  the  acoustic  pressure  (hammer  pressure)  for  the  vertical 
speeds  studied.  On  the  other  hand  current  structural  design  criteria  for  high-speed  craft  treat 
slamming  as  static  uniformly  distributed  pressures  with  considerably  lower  pressure  magnitudes 
(e.g.,  DNV  [2],  ABS  [3],  Lloyds  f4]).  This  raises  the  question  whether  the  structures  designed  and 
manufactured  by  those  criteria  are  conservative,  over  or  under  designed,  or  just  plainly  incorrectly 
designed.  The  present  study  tries  to  shed  some  light  on  this  complex  problem  by  analytically 
studying  a  simplified  model  of  a  bottom  panel  subjected  to  a  non-uniform  pressure  distribution 
traveling  at  various  speeds.  More  advanced  studies,  as  well  as  correlation  with  experimentally 
measured  slamming  response  using  the  Numerelte  research  craft,  are  underway. 

Early  analytic  research  on  slamming  was  done  by  von  Karman  [5]  using  a  momentum  approach, 
and  by  Wagner  [6]  using  two-dimensional  non-viscous  incompressible  flow.  Cointe  and  Armand 
[7]  studied  the  problem  of  an  impacting  cylinder  and  considered  nonlinearity  of  the  local  jet  flow. 
Zhao  and  Faltinsen  [8]  and  Faltinsen  [9]  improved  the  solution  of  Wagner  using  a  boundary 
element  method  and  indicated  a  superposition  of  asymptotic  expansions  of  high  pressure  at  the 
spray  root  and  a  following  lower  pressure  distribution.  Faltinsen  [9,  13]  reported  that  hydroelastic 
effects  are  mainly  relevant  for  local  impacts  when  the  deadrise  angle  is  small  and  the  duration  of 
the  impact  is  shorter  or  comparable  to  the  structure's  natural  period.  A  conformal  mapping 
technique  was  used  by  Mei  et  al.  [10]  to  study  the  impact  pressure  on  a  two-dimensional  body. 
Wet  deck  slamming  was  studied  theoretically  by  Faltinsen  [11]  using  a  hydroelastic  beam  model. 
An  initial  structural  inertia  phase  and  a  subsequent  free  vibration  phase  were  identified.  An 
asymptotic  theory  showed  that  the  maximum  bending  stresses  are  proportional  to  an  effective 
drop  velocity  and  are  not  sensitive  to  the  curvature  of  the  wave  surface  or  where  the  waves  hit  the 
beam. 
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Simply  stated,  slamming  appears  to  consist  of  an  initial  slamming  load  arriving  phase  followed  by 
a  vibration  phase.  In  this  paper  the  bottom  response  during  the  slamming  load  arriving  phase  is 
analytically  studied.  The  boat  bottom  is  modeled  as  a  one-dimensional  linear  elastic  Euler- 
Bemoulli  beam.  The  slamming  pressure  is  modeled  as  a  high  intensity  peak  followed  by  a  lower 
constant  pressure,  traveling  at  constant  speed  along  the  beam. 

Fluid-structure  interactions  are  at  this  time  ignored,  or  simply  included  as  a  constant  added-mass 
term.  Air  entrapment  which  may  have  a  large  influence  at  lower  deadrise  angles  is  also  ignored. 
The  deformation  is  assumed  to  be  sufficiently  small  that  linear-elastic  beam  theory  is  valid,  and 
that  the  geometry  of  the  deformed  bottom  is  not  significantly  different  from  the  undeformed  one. 
The  assumption  that  the  load  travels  at  constant  speed  across  the  beam  in  essence  implies  that  the 
vertical  velocity  of  the  boat  bottom  is  constant  during  the  slamming  event,  which  depending  on 
boat  particulars  (geometry,  mass,  etc)  may  or  may  not  be  a  reasonable  assumption.  The  equations 
are  solved  by  using  a  Fourier  sine  integral  transformation  in  space  and  a  Laplace-Carson  integral 
transformation  in  time,  as  done  by  Fryba  [12].  The  structural  response  during  the  slamming  load 
arriving  phase  is  given.  The  effect  of  slamming  load  traveling  speed  on  structural  responses  is 
presented. 


2.  Simplified  Analytical  Model  of  Slamming,  Two-Step  Load  on  Beam 

Consider  the  system  shown  in  Figure  1  in  which  a  bottom  panel  is  modeled  as  a  simply  supported 
Euler-Bemoulli  beam  subjected  to  the  slamming  load  q(x,t)  which  moves  with  constant  velocity  c 
from  one  end  to  the  other.  The  model  may  be  reasonably  realistic  for  bottom  panels  which  are 
long  relative  to  their  width,  as  is  common  in  boat  bottoms.  Simply  supported  edges  were  chosen 
as  a  reasonable  approximation  of  the  bottom  panels  of  the  Numerette  research  craft.  Its  bottom 
consists  of  ten  sandwich  panels  whose  cores  taper  off  and  vanish  by  the  edges,  resulting  in  a 
single  skin  "collar"  along  the  perimeter  of  each  panel.  The  thin  single  skin  collar  is  considerably 
more  compliant  in  bending  than  the  thick  sandwich,  thus  modeling  the  edges  as  simply  supported 
is  presumably  a  decent  approximation.  It  may  or  may  not  be  a  feasible  approximation  for  the 
edges  of  bottom  panels  in  other  boats.  The  deflection  of  the  beam  is  w(x,t),  where  x(0<x<Z)is 
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the  position  within  the  beam  and  t  (0  <t  <  L/c)  is  time.  In  Figure  1,  L  is  the  length  of  the  beam 
and  is  the  length  of  the  high-intensity  pressure  peak.  The  analysis  in  this  paper  terminates  at 
t=Llc  when  the  pressure  peak  reaches  the  right  end  of  the  beam  (corresponding  to  the  slamming 
pressure  reaching  the  chine  of  the  boat). 


Figure  1.  Left:  cross-section  of  boat  during  slamming,  showing  the  moving  slamming  pressure  q(x,t).  Right: 
Simplified  model  where  the  bottom  panel  is  represented  by  a  beam  and  the  slamming  pressure  q(x,t)  as  two  constant 

pressures,  q\  and  q2,  traveling  at  a  constant  speed  c. 


Using  the  Euler-Bernoulli  beam  assumptions  the  governing  equation  is 


d4w(x,t) 

dx4 


* 

+  M 


d2\v(x,t) 

~d? 


=  q(x,t) 


(1) 


where  El  is  bending  stiffness  (assumed  constant),  q(x,t)  is  load  per  unit  length,  p  is  total  mass 

(mass  of  beam  plus  some  added  mass  of  water)  per  unit  length  of  the  beam  (also  assumed 
constant).  The  added  mass  of  a  submerged  bottom  panel  is  usually  assumed  to  correspond  to  the 
mass  of  a  half  cylinder  of  water  with  diameter  L  and  length  d  (in  the  present  case  d  is  the  width  of 
the  beam).  The  total  added  mass  is  thus  pvnl}dl  8.  If  it  is  assumed  that  this  mass  is  evenly 
distributed  along  the  length  of  the  present  beam,  then  the  added  mass  per  unit  length  becomes 
pwnLd  1 8.  This  would  be  for  a  fully  submerged  panel,  whereas  at  the  beginning  of  the  slamming 
event  the  bottom  panel  is  essentially  dry  and  there  is  no  added  mass  term.  For  this  reason  a 
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constant  k  e  [0,l]  was  introduced,  such  that  ju*  =  jj.  +  kpwjiLd  / 8 .  For  a  fully  submerged  bottom 

panel  k=  1,  whereas  for  a  dry  panel  k= 0;  however  in  the  analysis  below  k  is  assumed  constant 
during  the  whole  slamming  event  (from  the  time  the  keel  touches  the  water  surface  to  the  time 
when  the  slamming  pressure  reaches  the  chine).  It  may  be  plausible  to  believe  that  the  two  cases 
k=0  and  1  in  some  sense  bound  the  behavior  of  a  bottom  panel. 


In  the  present  study  the  slamming  load  is  simplified  as  a  moving  step  load  of  the  following  form: 


q(x,t)  =  ql[\-H(x-ct)]-(ql-q2)[\-H(x-c(t-ll/  c))]  (2) 


where  q\  is  the  load  per  unit  length  of  the  initial  load  peak,  cj2  is  the  load  per  unit  length  of  the 
residual  load  following  the  peak,  and  H(x)  is  the  Heaviside  step  function, 


H(x)  = 


for  x<0 
for  x  >  0 


(3) 


Since  equation  (1)  is  linear,  superposition  applies  and  the  problem  can  be  solved  in  two  parts.  The 
first  set  of  equations  is 


_r34w(x,/)  *  d2w(x,t) 

El - +  p 

3r4  dt 2 


qi[\-H(x-ct)] 


(4) 


The  corresponding  boundary  conditions  and  the  initial  conditions  for  equation  (4)  are 


w(0,  /)  =  0  w(L,  t)  =  0 


(5) 
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d2w(x,  t) 


dx2 


=  0 


x=0 


d2w(x,t ) 


8x 2 


=  0 


x=L 


and 


w(x,0)  =  0 


dw(x,  t ) 


dt 


=  0 


t=0 


(6) 


The  second  set  of  equations  is 


EI^Sr^+v  =  -to  -  ?2)[i  -  H{x  -  c(t  -  /,  /c»] 

or  dt 


(7) 


The  corresponding  boundary  conditions  and  the  initial  conditions  for  equation  (7)  are 


w(0,t)  =  0  w(L,  t)  =  0 


(8) 


d2w{x,  t) 
dx2 


x=0 


=  0 


d2w{x,t ) 
dx2 


x=L 


=  0 


and 


w(^OUl/c=o 


dw(x,t ) 


dt 


=  0 


(9) 


If  w(x,f)  is  the  solution  to  the  slamming  problem,  wx (x, t)  is  the  solution  to  eqs.  (4-6)  and  w2(x,t) 
is  the  solution  to  eqs.  (7-9),  then 
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w(x,  t )  =  (x,  t )  +  w2  (x,  t ) 


(10) 


3.  Nondimensionalization  of  the  Problem 

In  order  to  reduce  the  number  of  parameters  the  following  dimensionless  quantities  are  introduced 


-  x 
x  =  — 
L 


r=i 

1  L 


—  w 
w  =  — 
L 


-  I) 

q  =  —  q 
El 


t=t 


I  El 

m'l4 


Ui  L2  - L 

c  =  c^—£f  M{x,t)=—M(x,t) 


(11) 


The  dimensionless  versions  of  eqs.  (4-6)  are 


dx  dt 

w(  0,  t)  =  0 

d2w(x,  t ) 


dx 2 


=  0 


q^-H(x 

-co] 

(12) 

w(l  ,t)  =  0 

d2w(x,  t ) 

_ o 

=  0 

(13) 

*=o 


dx 


x=l 


and 


w(x,0)  =  0 


dw(x,  t) 


dt 


=  0 


1=0 


(14) 


while  the  dimensionless  versions  of  eqs.  (7-9)  are 


+£2&a__(- 

dx  dt 


(15) 
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w(0,t)  =  0 


w(\,t)  =  0 


(16) 


d 2  w(x,  t ) 


dx 


=  0 


x=0 


and 


w(x,t) 


_  =  0 

t=li!c 


d2w(x,  t ) 


dx 2 


=  0 


X=1 


dw(x,  t ) 


dt 


=  0 


(17) 


t=h!c 


4.  Solutions  to  the  Equations 

For  equations  (12-14),  the  solution  wi(x,/)  can  be  obtained  by  using  a  Fourier  sine  integral 
transformation  in  space  and  a  Laplace-Carson  integral  transformation  in  time  [12].  Multiplying 
eq.  (12)  by  sin  jnx ,  where  j  is  an  integer,  and  integrating  with  respect  to  x  between  0  and  1, 
using  the  boundary  conditions  (13)  and  relations  (A1-A3)  in  Appendix  A,  the  following  is 
obtained: 


/ WC/,f)  +  W(j~t)  =  -^-(1  -  cos  jnct)  j  =  1,2,3  -  (18) 

J71 


where  W(j,t )  is  the  Fourier  sine  integral  transform  of  the  original  W\(x,t) , 


i  _ 

W  ( j ,  /)  =  |  w\  (x,  t)  sin  jnxd.  x  j  =  1,2,3  •  •  •  (19) 

o 
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(20) 


w\  (x,  t )  =  2^  W  (j,  t)  sin  j'kx 

7=1 

Set  a  =  j2n 2  and  6  =  jnc ,  then  equation  (18)  can  be  reorganized  as  follows 


W  (j~t) +a2W  (J~t)  =  cos  bt)  j  =  1,2,3-”  (21) 

J* 


Next  we  apply  the  Laplace-Carson  integral  transformation  on  equation  (21),  i.e.  multiply  each 

term  in  eq.  (21)  by  e~st ,  integrate  with  respect  to  t  between  0  and  oo ,  and  multiply  by  s  which  is 
a  variable  in  the  complex  plane.  Using  the  initial  conditions  (14)  and  relations  (B1-B3)  in 
Appendix  B,  we  can  get 


W\j,s) 


£i  bj_ 

jn  ( s 2  +  a2)(s2  +b 2) 


(22) 


where  W\j,s)  is  the  Laplace-Carson  transform  of  the  original  W(j,t), 


W'U,s)  =  s\W(j,~tyf,dt  (23) 

0 

Applying  the  inverse  Laplace-Carson  transformation  on  equation  (22),  using  relations  (B2)  (B4), 
the  solution  is  obtained, 


wu,~t)= 


<h 


— - j(b2  cos  at -a1  cos  bt) 

a  —b 


(24) 
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This  solution  is  not  valid  if  a=b,  i.e.  when  j  =  7ic .  However,  this  can  be  avoided  by  a  slight 
change  in  c ,  and  therefore  will  be  tacitly  ignored. 

With  equation  (20)  and  (24),  we  shall  get  the  solution  to  eqs.  (12-14),  where  0  <  /  <  1/c  and 
0<x<l, 


-  -  2 q.  A  1 

wi(x,t)  =  — fZ~ 


n 


y=i  y 


1+'  , 1  ,.Jt 

-\ 

r  cos  at  -  a  cos btj 

a2  -b2  ' 

sin  jnx ,  a  =  fn2,  b  =  jnc 


(25) 


The  solution  to  the  second  set  of  equations,  eqs.  (15-17),  is  obtained  by  a  time  shift  and  scaling, 


wi(x,t)  = 


0 


2(?2  1 


n 


24 


J= i  J 
2 


1  + 


a2-b2 


(b2  cos  a(t 


-a2  cos  b(t-=) 
c 


sin  jnx 


(oslsi) 

c 


C  C 


(26) 


Thus,  the  final  solution  for  a  moving  step  load  is 


w(x,t)  = 


2qx 


24 


^  M  f 

Av_L 

5  2-a  -5 

n  mJ 


1  +  ..  -  (b2  cos  at -a2  cosZtf) 


sin  jnx 


<h+  2  ,2 

a  —  b 


—  2  (b2  cos  at -a2  cos  bt) 


a2 -b2 


b2  cos a(t-=)-a2  cos b{t-^j 
c  c  , 


sin  y;rx 


(0<r<i) 


(i<t<i) 

c  c 


(27) 


With  the  deflection  known  other  quantities  of  interest,  like  the  bending  moment  or  the  shear 

forces  in  the  beam,  are  easily  obtained.  In  particular,  the  bending  moment  is 
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M  =  -El 


d2\v(x,t) 
dx 2 


(28) 


and  thus 


,  1  ,  (i b 2  cos  at  -a2  cos&n 
a2-b2V  ’ 

b2  cos  a{t  — hr)-a2  cosb{t-=)  I  sin  jnx 


{b2  cosat-a2  cosbt)  sin  jnx 


cos  at -a 


cos  at -a 


(o  <(<£-) 


c 


(29) 


— — b2  cos  a(t 


c 


c 


c  c 


In  the  next  section  a  slightly  different  pressure  profile  will  be  discussed. 


5.  Alternative  Pressure  Distribution,  the  Point-Step  Load 

As  mentioned  measurements  of  bottom  pressures  in  high-speed  boats  indicate  that  there  is  a  high 
pressure  peak  spreading  rapidly  over  the  bottom,  followed  by  a  considerably  lower  pressure  over 
a  large  area  of  the  bottom.  The  pressure  peak  was  represented  by  the  pressure  (times  width)  q\ 
above.  The  importance  of  the  pressure  peak  can  be  further  studied  by  representing  it  by  a  moving 
point-load,  F&,  followed  by  a  constant  pressure.  This  will  be  called  the  "Point-Step  Load"  in  the 

remainder  of  the  paper.  The  normalized  slamming  load  q(x,t )  can  then  be  expressed  as 


(30) 


where  §(x)  is  a  Dirac  pulse  and  qr  is  a  residual  pressure  following  the  point  load.  In  order  to 
compare  this  load  with  the  two-step  load,  the  point-load  load  is  set  equal  to  the  total  load  of  the 
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initial  peak,  and  the  total  load  when  the  slamming  reaches  the  right  end  of  the  beam  is  made  the 
same  for  the  two  different  loads;  thus 


Fo=qxh  qr=q2^-h) 


(31) 


The  governing  equation  for  this  load  is 


-  +  ~  =  Fd S(x -ct)+qr[\-H(x- ct )] 

dx  dt 


(32) 


together  with  the  boundary  and  initial  conditions,  eqs.  (13-14).  The  problem  is  solved  in  a  similar 
fashion  as  for  the  two-step  load,  and  with  the  help  of  (A4),  (B4)  and  (B5)  the  solution  becomes: 


Wi  (x,  t )  =  In  Fd 


1  . 


1  . 


P  a2  - b 2 
+  hi_f±- 

x5  jjj*l  a2  -  b2 


sin^t - sin  at 

b  a 


sin  jnx 


1-1 — — -( b 2  cos  at  -a2  cos  bt) 


sin  jnx 


(33) 


Mi  (x,  t )  =  2  n%  Fd 


J 


1  . 


1  . 


r  —sin  bt - sin  at 

pi  a  -b  \b  a  j 

i  r  i  / 


sin  jnx 


2 q  \  1  1  tj  2  —  2  i 

+  — 1+— 2 - cos  at-a  cos  bt) 


M  P Fl  a  —b 


sin  jnx 


(34) 


With  the  definition  of  the  point  load  and  the  residual  pressure  of  eq.  (3 1 ),  the  total  load  at  the  time 
the  slamming  reaches  the  right  end  of  the  beam  is  the  same  for  this  load  as  for  the  two-step  load. 
However,  at  any  other  time  the  total  force  is  higher  for  the  present  pressure  distribution. 
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6  Convergence  and  Error  Analysis 


In  this  section  the  errors  introduced  by  terminating  the  infinite  series  solutions  of  eqs.  (27)  and 
(29)  are  studied.  With  some  rearrangement  of  eq.  (25)  the  following  is  obtained, 


( V 

< 

2q] 

OO  I 

V  1  4- 

2qxc 

00 

V 

1 

/  \ 

1 

1  c? 

CM 

00 

v 

1 

n 5 

lu  ,5  + 

7-1  J 

7T5 

7=1 

.5/ .2  2  _2 

7  \7  ^  -c 

T 

n 3 

2-i 

7=1 

< 

2q] 

A- 1  1 

V  1  + 

2qxc 

Ji-1 

v 

1 

/  \ 

2^1 

j,-i 

v 

1 

n5 

Zj  {5  + 
7=1  J 

7T5 

7=1 

■5I  -2  2  — 2 

J  V  ^  -c  J 

l 

;r3 

2-r 

7=1 

/(/•v  -?) 

+ 


2q 


i+£lvj_ 


n 


Z4+ 


J-Ji  J 


2  qx 


71 


1 

Zi 

j=A  J 


(35) 


where  J\  is  an  integer  equal  or  greater  than  Vl  +  c2 1  tv  (so  j2n2  -c  >1).  Since  a  ^-series  with 

p>  1  is  convergent,  the  solution  w\(x,t)  is  also  convergent  (except  if  j  =  nc,  which  as  mentioned 
previously  is  disregarded). 

For  numerical  calculations  a  truncated  series  is  used, 


W 1 


1+— r- — x(b2  cos  at -a2  cos  bt) 
~J"  L  a2  -b  V  { 


Sm  J7TX 


(36) 


with  an  error  limited  by 


\error\ .  = 

I  \d 


2<?i  ^  1 


2ff](l  +  c  )  ^ J_  + 


1+— — T(b2  cosat  -  a2  cosbt) 
a  -b  ’ 


sin  jtux 


71 


i=J,  J 


2  qx 


n 


00  1 

J~J\J 


(37) 
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Since 


(38) 


the  error  estimate  becomes 


\error\ . 

I  1  a 


2n5{Jx  -l)4 


^3(J,-1)2 


(39) 


Similarly,  for  the  bending  moment  M i  =  - 


d2w\  (x,t) 


dx 2 


,  we  can  get  the  following  result 


*  7^9 


1h — — — ~(b2  cos  at -a2  cosbt] 
a  -b  ' 


sin  jtcx 


(40) 


where  Ji  is  an  integer  equal  or  greater  than 


l  +  Vl  +  4;rV 


-2 


I  In2 .  An  upper  bound  of  the  error  of 


the  bending  moment  M\ (x,t)  is 


<h 

1 

2  qxc 

1 

2qx 

c/  f  Ur  S 

1 

-l)2 

T 

3 a3(J2  -l)3 

(41) 
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7  Range  of  Parameters  for  Slamming 


Slamming  calculations  will  in  this  paper  be  performed  for  parameters  that  are  relevant  for  boats 
and  ships.  The  speed  that  the  slamming  load  travels  over  the  bottom,  c,  can  be  estimated  by 


Vn  1  Vn 

c  = - = - 

2tan p  cos/3  2sin {3 


(42) 


where  V  is  the  vertical  velocity  of  a  2D  wedge  dropped  in  water  and  (3  is  the  deadrise  angle  of 
the  bottom,  e.g.,  Faltinsen  [13].  It  will  be  assumed  that  deadrise  angles  of  the  boats  of  interest  are 
in  the  range  5°  to  45° ,  and  the  vertical  velocity  in  the  range  1  m/s  to  10  m/s.  Drop  tests  with  these 
parameters  have  been  performed  in  [1],  [14]  and  [16].  Using  these  values  the  speed  that  the 
slamming  load  travels  over  the  bottom,  c,  is  estimated  by  eq.  [42]  to  vary  from  2  m/s  to  200  m/s. 
Peak  pressure  and  duration  were  deduced  from  among  other  sources  the  drop  tests  of  [1],  [14]  and 
[16].  The  peak  pressure  range  and  peak  duration  in  Table  1  appears  to  cover  the  majority  of  such 
tests.  The  ratio  of  pi/p2  is  naturally  not  clearly  defined  from  experiments  since  real  slamming 
loads  differs  from  the  two-step  load  presently  assumed.  However,  for  the  purpose  of  presenting 
results  the  pressure  ratio  p\!pi  was  assumed  to  be  in  the  range  of  2~20. 


Table  1.  Range  of  parameters  studied 


Parameters 

Traveling 
load  speed  c 

[m/s] 

Peak  pressure 

[kPa] 

ratio  of 

P\!pi 

Duration  of 
slamming  load  peak 

[ms] 

Value  Ranges 

2-200 

10-8,000 

2-20 

0.01-2 
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Regarding  bottom  stiffness,  the  range  can  be  estimated  to  vary  from  that  of  a  very  soft  bottom 
panel  such  as  a  one  meter  wide  3  mm  thick  aluminum  plate  (E=7 0  GPa,  density  2700  kg/m3),  to  a 
stiff  bottom  such  as  a  one  meter  wide  sandwich  panel  with  two  1 5  mm  thick  carbon  fiber  skins 
(£■=100  GPa,  density  1500  kg/m3)  on  each  side  of  a  70  mm  thick  high  density  foam  core 
(negligible  stiffness,  density  250  kg/m3).  This  results  in  a  bending  stiffness  ( El)  ratio  of  the  soft 
bottom  to  the  stiff  bottom  of  approximately  3  x  10~5 .  With  this  ratio,  the  ranges  of  interest  of  the 
dimensionless  parameters  given  in  Table  2  are  obtained  (from  eq.  (11)).  The  length  1\  was 
estimated  as  traveling  load  speed  times  duration  of  the  pressure  peak,  resulting  in  the  normalized 

length  l\  of  Table  2  if  the  length  of  the  beam  is  1  m. 


Table  2.  Range  of  interest  of  the  dimensionless  parameters 


Parameters 

Traveling 

Peak  pressure 

Ratio  of 

Peak  load  length 

load  speed 

_ 

_  _ 

c 

9x 

h 

Value  Ranges 

0.01-320 

o 

d, 

X 

o 

2~20 

2x10_5~0.4 

In  the  next  section  dynamic  deflection  w  and  bending  moment  M  will  be  presented.  They  will 
be  normalized  by  the  maximum  static  deflection  wms  and  the  maximum  static  bending  moment 
M ms  that  result  if  the  beam  is  statically  subjected  to  a  two-step  load  with  the  same  properties  as 
the  dynamic  load  (i.e.,  as  in  Figure  1).  Let  x  =  a  be  the  right  edge  of  the  pressure  pulse  qx .  Then 


Wms  = 


max_ 

0<o<l,0<Jf<l 


Mms  =  max_ 

0<a<l,0<x<l 


M. 


(43) 
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where  ws  and  M s  are  the  normalized  static  deflection  and  bending  moment,  respectively;  please 
see  Appendix  C  for  more  details. 


8  Results  and  Discussion 

A  sample  calculation  will  reveal  some  of  the  features  of  dynamically  loaded  bottom  panels. 

Consider  a  simply  supported  beam  subjected  to  a  moving  slamming  load  with  the  following 

----  .  .  -1131 

parameters:  c  =  5 ,#,/ q2  =  5,  1 1  =0.01.  The  deflections  at  four  instances  (/  =  —=•,—=•,—=•  and  — ) 

4c  2c  4c  c 

are  shown  in  Figure  2,  and  the  bending  moments  are  shown  in  Figure  3.  The  vertical  axes 

represent  the  deflection  ratio  w/  wms  and  the  bending  moment  ratio  M  /  M  ms ,  respectively.  In  this 
example  the  Point-Step  Load  predicts  slightly  higher  deflection  and  bending  moment  of  the  beam 
in  most  instances,  as  would  be  expected.  In  this  case  the  maximum  deflection  reaches 
approximately  55%  of  the  maximum  static  deflection.  The  fact  that  the  load  is  moving  thus 
reduced  the  maximum  deflection  in  this  case.  This  is  not  generally  true  as  will  be  seen  shortly. 
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Figure  2.  Deflection  ratio  w/ w„,s  of  the  beam  under  moving  slamming  load 
with  c  =  5  ,  <7, 1  q2  =  5,  l\  =  0.01 
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Figure  3.  Bending  moment  ratio  M  /  M ms  of  the  beam  under  moving  slamming  load 


with  c  =  5  ,  qx  I  q2  —  5 ,  l\  =  0.01 


The  effects  of  the  slamming  load  traveling  speed  on  structural  responses  are  sketched  in  Figures  4 

and  5  when  l\  =0.01.  They  show  the  maximum  absolute  value  of  the  dynamic  deflection  and 
bending  moment. 


Wmd  =  _m_ax_ 

W 

Mmd  =  .max. 

M 

0</<1/c,0<jc<1 

0<r<l/c,0<jf<l 

(44) 


normalized  by  the  maximum  static  deflection,  Wm  ,  and  the  maximum  static  bending  moment, 

Mms ,  respectively.  Note  that  wmd  and  M„,d  depend  on  c ,  while  wms  and  Mms  do  not.  The 
results  in  Figures  4  and  5  demonstrate  an  important  phenomenon  of  the  structural  response  during 

the  slamming  load  arriving  phase.  When  the  dimensionless  slamming  load  moving  speed  is 
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relatively  low,  the  maximum  dynamic  deflection  and  moment  are  close  to  their  static 


counterparts,  as  expected.  However,  for  the  present  case  with  l\  =  0.01 ,  when  c  increases  to 
around  2  the  maximum  dynamic  deflection  is  approximately  50%  higher  than  the  maximum  static 
deflection.  The  same  is  true  for  the  bending  moment,  the  dynamic  bending  moment  is 
approximately  50%  higher  than  the  static  one.  There  is  a  form  of  resonance  occurring  in  the 

structure.  Considerable  vibrations  are  occurring  during  the  slamming  load  arriving  phase  when  c 
is  relatively  low.  From  equations  (1)  and  (8)  the  eigenfrequencies  of  a  simply  supported  Euler- 


Bemoulli  beam  can  be  easily  found:  cot  = 


2-2 
K  X 


(i  =  1,2,3  •  •  •) .  If  we  define  a  characteristic 


velocity  as  c,  = - then  the  dimensionless  characteristic  velocity  will  be  c,  = - .  The  first 

2  n  2 

three  dimensionless  characteristic  velocities  are  1.6,  6.3  and  14.1.  The  two  lower  graphs  in 
Figures  4  and  5  indicate  that  when  l\  =  0.01  the  most  severe  structural  response  occurs  when  the 

propagation  speed  c  is  slightly  higher  than  the  first  characteristic  velocity.  Faltinsen  [11]  also 
pointed  out  that  maximum  strains  occur  during  the  free  vibration  phase  and  mainly  the  lowest 
eigenmode  is  of  importance  at  the  time  scale  when  maximum  strains  occur. 


Presumably  due  to  inertia  effects,  at  higher  propagation  speed  c  the  maximum  deflection  and 
bending  moment,  wmd  and  M„,d ,  decrease  rapidly.  When  c  is  larger  than  the  third  characteristic 

velocity,  14.1,  the  maximum  structure  response  ratios,  wmd  /  wms  and  MmdIMms,  are  rather 
small.  This  is  referred  to  as  the  inertia  phase  by  Faltinsen  [11],  and  implies  that  the  slamming 
force  is  essentially  balanced  by  structural  inertia  forces.  After  the  inertia  phase,  the  structure  starts 

to  vibrate  with  an  initial  velocity  obtained  at  the  end  of  the  inertia  phase.  For  c  larger  than  40,  the 
maximum  dynamic  deflection  and  bending  moment  at  the  end  of  the  slamming  load  arriving 
phase  is  less  than  10%  of  the  corresponding  maximum  static  values. 


Some  effect  of  varying  the  peak  pressure  ratio,  qx!  q2,  is  also  illustrated  in  Figures  4  and  5.  With 
h  =  0.01 ,  increasing  the  pressure  ratio  from  1  to  20  leads  to  higher  deflection  and  bending 
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moment  ratios,  wma  lwms  and  Mma  /Mms.  In  other  words,  a  slamming  load  with  the  same  peak 
pressure  but  higher  pressure  ratio  will  result  in  a  more  severe  structural  response  compared  with 
the  corresponding  static  response. 


Figure  4.  Maximum  deflection  ratio  versus 
speed  of  the  slamming  load  for  five  different 

pressure  ratios  qx  /  q2  when  l\  =  0.01 .  The 

lower  graph  is  a  zoomed-in  version  of  the 
upper  graph. 


Figure  5.  Maximum  bending  moment  ratio 
versus  speed  of  the  slamming  load  for  five 

different  pressure  ratios  qx/q2  when 

/ i  =  0.01 .  The  lower  graph  is  a  zoomed-in 
version  of  the  upper  graph. 
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An  effect  of  peak  load  length,  l\ ,  is  shown  in  Figures  6  and  7  for  ql  /  q2  =  5.  Three  different  h 

were  used,  0.001,  0.01  and  0.1.  The  deflections  and  bending  moments  were  normalized  by  their 
static  equivalents  as  before  in  this  paper.  The  figures  indicate  that  both  the  normalized  maximum 

deflection  and  dynamic  bending  moment  increase  with  increasing  peak  load  length,  l\ .  Hence, 
slamming  load  with  long  duration  of  the  peak  pressure  appears  to  generate  more  deflection  and 
bending  moment  compared  with  the  static  equivalents. 


c 

Figure  6.  Maximum  deflection  ratio  versus 
speed  of  the  slamming  load  for  three  different 


l\,  when  ql/q2=5 


c 

Figure  7.  Maximum  bending  moment  ratio 
versus  speed  of  the  slamming  load  for  three 


different  l\ ,  when  qi  /  q  2  =5 


Some  further  insight  into  slamming  can  be  gained  by  studying  the  influence  of  the  total  force  of 
the  initial  peak  pressure.  The  Point-Step  load  outlined  previously  was  used  for  this  purpose,  i.e.,  a 
point  load  Fd  preceding  a  constant  pressure  qr.  The  maximum  dynamic  deflection  and  the 
maximum  dynamic  bending  moment  were  calculated  for  different  point  loads  Fd-  In  this  case  the 
deflection  was  normalized  by  the  static  deflection  of  a  beam  loaded  by  a  distributed  pressure  qT 


only,  i.e.,  by  ws  =  / 3 84  .  Likewise,  the  dynamic  bending  moment  was  normalized  by 


Ms  = 


qr! 8  .  The  results  are  presented  in  Figures  8  and  9  where  deflection  and  bending  moment 


are  plotted  versus  FD  /  qrL . 
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Figure  8.  Maximum  deflection  versus  point 
force  ratio,  normalized  by  maximum 
deflection  from  static  evenly  distributed  load 
qr. 


Figure  9.  Maximum  bending  moment  versus 
point  force  ratio,  normalized  by  maximum 
moment  from  static  evenly  distributed  load  qr. 


Figure  8  and  9  show  that  the  maximum  deflection  and  bending  moment  of  the  structure  increase 
monotonously  with  increasing  point  force  ratio,  FD  /  qrL ,  as  expected.  Further,  as  previously  seen 

when  c  is  close  to  the  first  characteristic  velocity,  1 .6,  the  maximum  deflection  and  the  maximum 
bending  moment  of  the  structure  are  the  largest.  After  this  high  deflection  and  high  bending 

moment  phase,  ww  and  M„,d  decrease  with  increasing  speed  of  the  slamming  load.  Further,  the 
figures  indicate  that  reasonably  small  point  loads  have  only  little  effect  on  the  structural  response 
when  the  slamming  load  travels  slowly.  For  example,  the  increase  in  deflection  and  bending 

moment  is  just  a  few  percent  (relative  to  ws  and  Ms)  for  FD  / qrL  =  0.1  and  c  =  0.1  or  c  =  1 . 
However,  at  higher  speeds,  such  as  c  =  2  or  c  =  5,  the  increase  is  quite  considerable.  For  c  =  2 
the  increase  in  deflection  and  bending  moment  over  and  Ms  is  approximately  12  %  while  for 

c  =  5  the  deflection  increased  31%  and  the  bending  moment  38%  in  spite  of  the  fact  that  the  total 
load  is  only  1 0  %  higher  because  of  the  point  load. 


23  (31) 


9  Conclusions 


The  initial  phase  of  slamming  was  studied  using  a  simplified  beam  model  subjected  to  moving 
loads.  At  slow  speed  of  the  moving  load  the  maximum  dynamic  deflection  and  the  maximum 
dynamic  bending  moment  do  not  differ  much  from  their  static  counterparts,  but  when  the  speed 
approaches  a  characteristic  speed  of  the  panel  the  maximum  deflection  and  bending  moment 
increase  on  the  order  of  50  %.  At  speeds  above  a  few  times  the  characteristic  speed,  the  maximum 
deflection  and  bending  moment  decrease  below  the  static  values  (during  the  initial  phase  under 
study).  Both  keep  decreasing  monotonically  as  the  speed  increases.  This  suggests  that  bottom 
panels  could  be  tailored  to  avoid  large  deflections  and  bending  moments.  In  particular  panel 
stiffnesses  could  be  avoided  for  which  the  characteristic  speed  of  the  panel  is  near  the  propagation 
speed  of  the  slamming  load.  Granted,  this  may  be  difficult,  or  even  impossible,  to  achieve  for  any 
speed  and  any  sea  state  for  a  particular  boat.  Nevertheless,  it  may  be  possible  to  design  a  boat's 
bottom  such  that  the  effects  of  the  most  severe  condition  (speed,  sea  state)  are  reduced.  It  should 
further  be  noted  that  the  speed  and  the  pressure  of  the  moving  load  depend  on  deadrise  angle  of 
the  hull;  for  example  reducing  the  deadrise  angle  would  increase  the  speed  of  the  slamming  load 
(which  may  be  beneficial)  but  also  the  pressure  (which  would  not  be  beneficial).  More  analyzes 
using  more  refined  models  would  be  required  to  gain  a  better  understanding  of  the  potential  of 
tailoring  bottom  panels  for  slamming. 


The  leading  edge  of  the  slamming  pressure  is  typically  characterized  by  a  high  pressure  peak.  If 
the  force  in  this  peak  is  on  the  order  of  10%  of  the  total  force  on  the  bottom  panel,  then  at  slow 
slamming  load  propagation  speeds  (less  than,  say,  half  the  characteristic  speed  of  the  panel)  this 
pressure  peak  does  not  increase  maximum  deflection  or  maximum  bending  moment  significantly. 
However,  at  higher  speeds  (on  the  order  of  three  times  the  characteristic  speed)  the  maximum 
deflection  and  maximum  bending  moment  increase  30-40%  due  to  this  pressure  peak  (compared 
to  the  response  from  a  static  evenly  distributed  load  with  no  pressure  peak),  which  is  very 
significant. 
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Appendix  A 


Fourier  sine  integral  transformation 


Equation 
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Appendix  B 

Laplace-Carson  integral  transformation 


Equation 
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Appendix  C 


The  dynamic  moment  is  normalized  by  the  maximum  static  moment.  The  maximum  static 
moment  within  the  beam  depends  on  where  the  (slamming)  load  is  applied.  Introduce  the  distance 
a  such  that  the  right  end  of  the  load  q\  in  Figure  1  is  located  at  x=a.  In  the  dynamic  case  a=ct.  The 
maximum  static  bending  moment  cannot  occur  for  a<l\  since  additional  load  will  increase  the 
maximum  moment.  For  a>l\  the  bending  moment  in  a  beam  that  is  statically  loaded  as  in  Figure  1 
is 


M  = 


q2 x2  (/i2<7i  —  2a/,  qx  +  2 Llx  qx  -  a2q2  +  2a/,  q2  -  l2q2  +  2 Laq2  -  2 Llx  q2  )x 

- b -  ,  0  <  x  <  a  -  L 

2  2  L  1 

qxx2  [l2q{  -2allql  +2Laq{  -a2q2+lalxq2-l2q2)x  [a  -  lxf  [qx  -  q2) 

- 1 - ,  a-L  <x<a 

2  2L  2  1 

\l2qx  -  2a/,  qx  -  a2q2  +  2a/,  q2  -  l2q2  Jbc  [l2qx  -  2a/,  qx  -  a2q2  +  2a/,  q2  -  l2q2  J 


2  L 


a<x<L 


The  maximum  moment  can  occur  either  when  a=l\,  when  a=L,  or  when  l\<a<L.  In  either  case, 
since  M(0)=M(L)=0  and  M  and  its  slope  are  continuous,  the  maximum  moment  will  be  found 
where  dMdx=0.  The  three  different  cases  are: 

Case  1,  a=l\ 


Ml 


g.x2  [  (2Z/,  -l2\x 
2  2  L 

A q\x  ,  h *?i 
2L  2  ’ 


,  0  <  x  <  /, 

/,  <x<  L 


The  maximum  moment  for  this  case  is 


Ml 


/,2(2L-/,)2g, 

8  L2 
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Case  2,  l\<a<L 


The  maximum  will  be  in  the  range  0 <x<a  since  for  x>a  the  derivative  dM/dx  is  never  zero.  The 
maximum  moment  will  be  found  in  the  range  x<a-li  if 


k  <k  “ 


aqx 


—  L(]\  -  aq2  +  -Ja^qf  -  2aLq\  +  L2qx  —  a2qx  q2  +  2aLqx 


<h 


9i  ~<h 


and  otherwise  in  the  range  a-h<  x<a.  The  maximum  moment  is 


Mmax  = 


(  k  (A  +  )(</j  q2^+  o  q2  +  2a{lx  [cjx  <72 )  ^q2 )) 

$lFq2 


if  lx  <  lx 


( 2alxqx  -l\qx  +a1q2  —2  alxq2  +  lxq2  j2alx  qx  -l\qx  -4aLqx  +  4L2qx  +a2q2  -2alxq2  +  lxq2) 


8  L2qx 


if  l  > 


Case  3,  a-L 


The  moment  is 


M3  =< 


q2x  (ji  q  i  k  Qi  L  q2^)x 
- +- - — 


0  <  x  <  L—lx 


2  2  L 

qx  x  {lx  qx  —2Llx  qx  +2L  qx  —  L  q2  +2Llx  q2  —lx  q2  (l—Ix  )  [qx  —  q2 ) 


-+- 


2  L 


L—lx  <x<L 


The  maximum  moment  will  be  found  where  dM/dx=0,  which  is  in  the  range  x<L-l\  if 


k  < 


Vizi  ~i~q2  +1 
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and  otherwise  in  the  range  x>L-l\.  The  maximum  moment  is 


M3  =< 

max 

(1i(q  i  ~<h  )+L2q2)2 

,  if 

%L2q2 

((2iZ,  -Z,2X<Z,  -?2)+i2?J 

8  L2q, 

0  <  /j  < 


if  L>lt  > 


-JqJq~2+\ 
L 


V2i  iQi 


Total 

The  maximum  moment  is  the  maximum  of  Ml,  A/2,  and  M3, 
Mms  =  ma  x[Mlm ,  M2^ ,  M3max  ] 


+  1 
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